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The nonequilibrium critical dynamics of the Ising magnet on a fractal substrate, namely the Sier- 
pinski carpet with Hausdorff dimension du =1.7925, has been studied within the short-time regime 
by means of Monte Carlo simulations. The evolution of the physical observables was followed at 
criticality, after both annealing ordered spin configurations (ground state) and quenching disordered 
initial configurations (high temperature state), for three segmentation steps of the fractal. The topo- 
logical effects become evident from the emergence of a logarithmic periodic oscillation superimposed 
to a power law in the decay of the magnetization and its logarithmic derivative and also from the 
dependence of the critical exponents on the segmentation step. These oscillations are discussed in 
the framework of the discrete scale invariance of the substrate and carefully characterized in order 
to determine the critical temperature of the second-order phase transition and the critical exponents 
corresponding to the short-time regime. The exponent 6 of the initial increase in the magnetiza- 
tion was also obtained and the results suggest that it would be almost independent of the fractal 
dimension of the susbstrate, provided that dn is close enough to d = 2. 
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I. INTRODUCTION 



The study and characterization of continuous phase 
transitions occurring on fractal substrates have attracted 
much attention during the past decades in the research 
fields of statistical physics and materials science. The 
very well-known general belief in universality states that 
for a given symmetry of the order parameter and range 
of interactions, the critical behavior only depends on the 
dimensionality d, i.e. the influence of the underlying 
structure becomes negligible at the critical point when 
the correlation length is much larger than the cell spac- 
ing. In particular, Ising models embedded in different 
lattices in the absence of external magnetic fields, with 
short-range interactions and for a given d >1, exhibit 
a continuous phase transition characterized by identical 
critical exponents, with the critical temperature decreas- 
ing to zero at the lower critical dimension d=l. Instead of 
the replication of an elementary cell by translation, frac- 
tals are constructed by the iteration of a generating cell, 
consequently the topological details of the generating cell 
are present at any scale. In this way, the critical expo- 
nents also depend on other geometric and topological pa- 
rameters^uch as the ramification and lacunarity of the 
fractal 111 (211. It has been shown that a second-order phase 
transition at nonzero temperature occurs only if the frac- 
tal substrate has an infinite ramification or 
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On the other hand, the lacunarity, which was introduced 
by Mandelbrot in order to give a measure of the devi- 
ation from the translational invariance can be ob- 
tained considering the departure from the power law of 
the mass scaling definition of the fractal Hausdorff di- 
mension (d//), namely m(i?) oc R'^" (where m{R) is the 
mass inside a sphere of radius R centered inside a frac- 
tal structure). Direct quantitative studies of the Ising 
model on fractal substrates have shown that the values 
of the ratio of critical exponents "f/v depend upon the 
lacunarity and indicate that by increasing the lacunarity 
one has to introduce higher corrections in the finite-size 
scaling associated with the critical behavior of physical 
observables Most of the previous work on this topic 
was based on the same class of fractals, built in the same 
way as the standard Sierpinski carpet (SC(b,c)), which 
have an infinite ramification order. A hypercube in d 
dimensions is segmented into b'* subhypercubes and c'^ 
of them are then removed, this segmentation process is 
iterated on the remaining subhypercubes. The mathe- 
matical fractal is obtained after an infinite number of 
segmentation steps (k). However, for practical purposes, 
a fractal is constructed by applying a finite number of seg- 
mentation steps. So, the physical observables show not 
only finite-size effects caused by the size of the system, 
but also a lower cutoff due to the segmentation step. As 
a matter of fact, results obtained by performing Monte 
Carlo simulations depend on the number of segmentation 
steps, so it is interesting to know the way these observ- 
ables converge to the ones expected for the mathematical 
fractal. The key point here is that in order to obtain crit- 
ical exponents from simulations in equilibrium, one has 
to apply finite-size scaling analysis that is hindered by 
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topological scaling corrections. These corrections and the 
critical slowing down increase when the fractal dimen- 
sion decreases towards one. Furthermore, due to these 
effects the segmentation steps studied could not be large 
enough to truly observe the asymptotic behavior. So, 
under these circumstances it is not surprising that also 
the critical exponents could not be easily obtained^, 0. 
Furthermore, Pruessner et al ^ have questioned the va- 
lidity of finite-size scaling on fractal substrates arguing 
that each segmentation step represents a new thermody- 
namic system that cannot be treated as a scaled version 
of the previous one. 

During the last years, nonequilibrium critical dynam- 
ics has been developed and a dynamic scaling form was 
found that is already valid in the short-time regime 
0, [Tol |. In addition to the dynamic exponent z, the 
static exponents, originally defined in equilibrium, also 
appear in the short-time scaling form. This fact leads 
to new methods for the numerical measurement of all 
the static and dynamic critical exponents, as well as the 
critical temperature. Since the measurements now are 
carried out within the short-time regime, they do not 
suffer from critical slowing down. Because of the small 
nonequilibrium correlation length, it is also easy to over- 
come finite-size effects. Due to the above-mentioned ad- 
vantages the nonequilibrium critical dynamics within the 
short-time regime is a promissory tool for the study of 
continuous phase transitions on fractal substrates. Hith- 
erto, the topological effects on the short-time dynamic 
behavior and the possible corrections to scaling have not 
been systematically investigated. Recently, we have re- 
ported soft oscillatory deviations from the power-law de- 
cay of the magnetization in the short-time regime, llj for 
the Ising model on the SC(3,1) with = 1.8927. On 
the other hand, logarithmic periodic oscillatory devia- 
tions in the behavior of physical observables have been 
reported for several systems that present fractal charac- 
teristics embedded into the dynamics of the model and/or 
the substrate. Some examples are the wave propagation 
in fractal svstems Il2l. the Blume-Capel model on the 
Sierpinski gasket p^. the dynamics of biological systems 
such as the bronchial tree0|, proteins [l^, magnetic and 
resistive effects on a system of wires connected along the 
Sierpinski gasket . and random walk through fractal 
environments^^ . 

The purpose of this paper is to study the occurrence of 
topological effects on the scaling laws that describe the 
nonequilibrium short-time critical dynamics of a mag- 
netic system on a fractal substrate. For this purpose we 
have performed extensive numerical Monte Carlo simu- 
lations of the Ising model on the SC(4,2). Three differ- 
ent segmentation steps are investigated in order to show 
how the lower cutoff due to the finite-segmentation step 
affects both the critical exponents and critical tempera- 
ture. This procedure also allows us to study the conver- 
gence of the measured values to the values expected for 
the mathematical fractal. 

In Section II we introduce the model, we give the sim- 



ulation details and remark on the main features of the 
short-time dynamic scaling. In Section III we present, 
analyze and discuss the short-time dynamic behavior of 
systems starting from both disordered (III A) and or- 
dered initial states (III B). Critical exponents obtained 
under these two conditions are compared. Also, in sub- 
section III C we discuss the topological characteristics of 
the substrate and their influence on the dynamic critical 
behavior of the physical observables. In Section VI we 
present our conclusions. 

II. THE MODEL AND THE METHODS 
A. The Ising Model on Sierpinski Carpets 

Figurenshows a sketch of the fractal substrate used in 
the present work, namely the SC(4,2), for the segmenta- 
tion step k = 3. The Hausdorff dimension of the math- 
ematical fractal is given by dn = '°^iog4^ —1-7925. In 
the case of the fractal, the convergence of relevant prop- 
erties towards the thermodynamic limit occurs when the 
structure is constructed by the iterative process. Such 
convergence is reflected in the dependence of the mean 
number of nearest neighbors per site on the segmentation 
step, which converges to a constant value as k goes to in- 
finity. In order to study the critical behavior by means of 
Monte Carlo simulations it is necessary to simulate large 
enough segmentation steps to ensure that the asymptotic 
region is reached. For the SC(4,2) the deviation of the 
mean number of next-nearest neighbors from the value 
corresponding to the thermodynamic limit (mathemati- 
cal fractal), using periodic boundary conditions and de- 
termined by means of the transfer-matrix method, be- 
comes negligible for k > 5 17J. 

As mentioned, the studied fractal has an infinite ram- 
ification order, which implies that the Ising model would 
exhibit a continuous phase transition at finite tempera- 
ture. In order to perform the simulation the spins were 
placed at the center of the occupied subsquares and con- 
sequently, the number of spins increases as a power law 
of the lattice size with an exponent given by dn- The 
Hamiltonian of the system is then given by 

H = -jJ2s.,Sj, (1) 

where Si assumes the values ±1, the sum runs over all in- 
teracting first-neighbor pairs and the positive exchange 
coupling constant J corresponds to the ferromagnetic 
case. 

B. Short-Time Dynamic Behavior at Second-Order 
Phase Transitions 

In this section we give a brief summary on the 
short-time dynamic (STD) method for characteriza- 
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length critical exponents, respectively. The scaling func- 
tion behaves as F{x) ~ 1 for x ^ and F{x) ~ - 
for X ^ oo. Extrapolating the results to mo = 0, the 
exponent 9 of the initial increase in the magnetization, 
which is a new nontrivial critical exponent, can be ob- 
tained. We remark that A/^(t) and A{t) have to be mea- 
sured starting from fully uncorrelated configurations with 
strictly zero magnetization, i.e for T = oo. One should 
also expect that the generalized hyperscaling relationship 
given by d^ff — ^ = 7 would hold for the Ising model 
on a fractal substrate (see also equation ©). 

For a dynamic relaxation from a completely ordered 
state (mo = 1) corresponding to the ground state at 
T = 0, annealed to Tc, the magnetization, the loga- 
rithmic derivative of the magnetization with respect to 
the reduced temperature r = "^Z^" , and the second- 
order Binder cumulant should follow power laws in time, 
namely 



FIG. 1: Sketch of the Sierpinski carpet SC(4,2) iterated up to 
the A: = 3 segmentation step. Spins are placed at the center 
of filled squares. 



tion of continuous phase transitions. According to 
STD arguments,^, a magnetic system at high temper- 
ature and with a small magnetization mo ^1 suddenly 
quenched to the critical temperature T^, presents a uni- 
versal dynamical evolution, which sets right after a mi- 
croscopic time scale imic^ large enough in the microscopic 
sense but still very small in the macroscopic sense. If 
during this time interval the nonequilibrium correlation 
length remains smaller than the system size the short- 
time dynamics is free of finite-size effects. According 
to the short-time dynamic scaling^,^^ the magnetiza- 
tion, its second moment and the autocorrelation function 
should follow power-law scaling relations given by 
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respectively. Here (...) denotes the average over samples, 
N is the total number of spins in the sample, is the ex- 
ponent of the initial increase in the magnetization, and 
de// is the scaling dimension whose value may be dif- 
ferent than djy. Also, (3 and v are the standard- usually 
defined in equilibrium-order parameter and correlation 
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respectively. For T ^ but within the critical region, 
the power-law behavior is modified by a scaling func- 
tion, which for the magnetization is given by M* (t~ t). 
So, this fact can be used to determine the critical tem- 
perature from the localization of the optimal power-law 
behavior. 

We have found that the dynamic evolution of the mag- 
netization when the system is annealed from T = to a 
higher temperature T close to the critical point is more 
sensitive to this final temperature than when the system 
is quenched from T=oo. Furthermore, due to the large 
initial value of M and its slow decrease upon annealing, 
the statistical fluctuations are less prominent, and there- 
fore it would be expected that the topological effects may 
be easier to detect. 



C. Monte Carlo Simulations 

In order to study the effects of the topology of the 
fractal on the STD and the convergent behavior of the 
critical exponents with the segmentation step, we car- 
ried out Monte Carlo simulations for segmentation steps 
A: = 4, 5, and 6 (system size L = 256, 1024, and 4096, 
respectively) using periodic boundary conditions. Sim- 
ulations started either from the complete ordered state 
or from a disordered state with zero or a small magne- 
tization. In the last case the initial magnetization mo 
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was settled by flipping, in a random disordered config- 
uration, a definite number of spins placed at randomly 
chosen sites in order to get the desired value of mo. The 
time evolution of the system was updated by means of a 
Metropolis algorithm usin g th e Marsaglia-Zanan pseudo- 
random number generator |l9j. The time unit, defined as 
a Monte Carlo step (MCS), involves attempts to update 
N randomly selected spins. The time evolution was fol- 
lowed, depending on the initial state, from 10^ up to 2 x 
10^ MCS. The magnetization, the autocorrelation, and 
the second moment of the magnetization were averaged 
over a number Us of samples with the same initial condi- 
tions but using different configurations. 



III. RESULTS AND DISCUSSION 

A. Dynamic Evolution from the Ordered State 

Figure |2 shows the decay of the magnetization as a 
function of time obtained at different annealing tempera- 
tures for the segmentation step k=5. As can be observed, 
log-periodic oscillations are present in the time evolution 
of the magnetization over the full range of time and for 
all temperatures. We believe that the log-periodic struc- 
ture exhibited by the data is a typical feature of systems 
with Discrete Scale Invariance (DSI)^3|- Section 
III C we first show that the topology of the substrate (the 
Sierpinki carpet in our case) possesses this property and 
subsequently, we formulate arguments linking the spatial 
characteristic length of the fractal to the time character- 
istic period of the dynamic behavior of the Ising system. 
Now, in order to account for our numerical observations 
we assume that the power-law behavior expected for the 
time evolution of the magnetization (equation (S)) should 
be replaced by 

M{t) oc t-^ (^ + A cos [j^jT) log(<) + ) ' (8) 

where log(P) is a logarithmic period, while A and (j) 
are the amplitude and the phase of the oscillation, re- 
spectively. An similar expression to equation® has 
been used by other authors to describe the time depen- 
dence of systems that exhibit a DSI in their dynamic 
behavior, such as the energy release on the approach of 
impending rupture, earthquakes 12], and biological sys- 
tems (proteins) . In order to study the critical regime 
we settle down the temperature range that exhibits the 
smallest deviations from a power-law behavior, namely 
1.09 < Tc < 1.13. Subsequently, we determine the criti- 
cal temperature by finding the smallest deviation of the 
data from equation |(SJ). For fc = 5, this study yields 
Tc = 1.11(1), where the error bar was assessed by consid- 
ering the closest pair of temperatures that present notice- 
able but small deviations. The inset of Figure|21shows the 
oscillatory component of M{t) that is nicely fitted by the 
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FIG. 2: Log-log plots of magnetization versus time obtained 
for diflFerent temperatures close to criticality, starting from 
ordered initial conditions (mo = 1) for the segmentation 
step k=5. The critical temperature Tc{k = 5) = Lll(l) 
was determined by obtaining the best fit of M{t) to equa- 
tion The inset shows the oscillating component of M(t) 

{Y = mty-ct^^ versus t in log-scale for T=Tc, T=L10 
ct ~ 

and T=1.12, which we fit with a cosine function (solid line). 
More details in the text. 



cosine function for T = 1.11, and also shows the depar- 
tures from the cosine function for T = 1.10 and T = 1.12. 
Logarithmic periodic oscillations are also observed for the 
segmentation steps k = A and 6, and the procedure pre- 
viously described was also used to determine the corre- 
sponding critical temperatures, as reported in Table Q 
(2nd column). Figure |3 shows the magnetization decay 
at the critical temperature for the segmentation steps 4, 
5, and 6. The fits were carried out after disregarding 
an initial time interval tmin = 30MCS, and the obtained 
parameters are reported in Tabled For k — 4, in order 
to assure that the noncquilibrium correlation length re- 
mains smaller than the system size, the considered time 
interval was SO-IO^'MCS. As observed in the inset of Fig- 
ure 01 the logarithmic periodic oscillations of the magne- 
tization decay have similar shape for k=5 and 6, but they 
are slightly different for k= 4. The critical temperature 
Tc{k = 6)= 1.10(1) is in relatively good agreement with 
determinations performed by means of Monte Carlo sim- 
ulations in equilibrium and obtained by applying finite- 
size scaling, namely <1.1780 and Tc =1.077(3)^]. 
Furthermore, a careful inspection of the data reveals a 
systematic but convergent decrease in both Tc{k) and 
P/i'z (TableQl 3rd. column) when the segmentation step 
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TABLE I: Critical temperatures and parameters obtained by 
fitting the time dependence of the magnetization according 
to equation ratio of critical exponents P/vz, and loga- 
rithmic period (log(P)), amplitude (A) and phase (0) of the 
oscillation. Data obtained by starting the simulations from 
an ordered initial state and for segmentation steps k= 4, 5, 
and 6. 



k T,(fc) 



A 



log(-P) 



4 1.14(1) 0.0145(1) 0.0023(2) 2.11(7) tt 

5 1.11(1) 0.0116(1) 0.0029(1) 2.15(1) tt 

6 1.10(1) 0.0110(1) 0.0029(1) 2.16(1) tt 
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FIG. 3: Log-log plots of magnetization versus time obtained 
at the critical temperature starting from ordered initial condi- 
tions (mo = 1) for the segmentation step k=4, 5, and 6. The 
corresponding fits are also shown by mean full lines. The inset 

Ct ~ ^ 



shows the shape of the oscillation of M{t) {Y = Mill- 
et 

versus t in log-scale. 



is increased, suggesting that our results could be taken 
as upper bounds. The exponent (i/yz for k = 6 is no- 
tably smaller than both the accepted value for the Ising 
model in d = 2, given by (3/vz — 0.0577(3) 2Lj, and our 
estimation for the SC(3,1), given by P/vz = 0.03412(7) 
In this way at the critical point a larger deviation 
from the translational symmetry, i. e. an increase in the 
lacunarity, makes the decay of the long-range order slow 
down. 

The logarithmic derivative of the magnetization with 
respect to r is evaluated by taking the difference between 
the values of M{t) at two temperatures close to the crit- 
ical one (Tc ± 0.01). The results of this calculation are 
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FIG. 4: Log-log plot of the logarithmic derivative of magneti- 
zation versus time obtained at Tc by starting from the ordered 
initial conditions (mo = 1). Results obtained for three differ- 
ent segmentation steps as listed in the figure. The full lines 
correspond to fits obtained for t > 100 MCS, according to 
equation Q for = 5 and 6, and to equationj^jlfor = 4, 
respectively. See details in the text. 



shown in Figure 0] where the logarithmic periodic oscil- 
lations can be clearly observed for the segmentation step 
/c = 6. For k = 5 the oscillations are also suggested by 
the behavior obtained for t > lO'^MCS. In order to ob- 
tain the exponent l/i'z, we propose the same correction 
for the power law as that appUed to the magnetization 
and given in equation ©, namely 



Vr{t) OC t- 



l-f i? cos 



27r 



log(P) 



log(t) + $ 



(9) 



where B and $ are the amplitude and the phase of the 
oscillations, respectively. The fits were performed within 
the range 10^ to 2x lO^MCS giving a value for log(P) — 
2.0(2) and the exponents l/vz listed in Table ITU For k — 
4, the presence of oscillations cannot be detected within 
the time range considered, i.e. from 10^ to lO^'MCS. In 
order to give a crude estimation for the critical exponent, 
we fitted the data with equation ^ (see Table^ . Figure 
|31 shows the determination of the exponent ^^^^ (see Ta- 
blem from the time dependence of the Binder cumulant, 
according to equation (|7|) . The trend of the data for the 
exponents l/vz and namely a convergent decrease 

when k is increased, is consistent with the previous obser- 
vations reported for SC(3,1)|^^3 and strongly suggests 
that our results should be taken as reliable upper bounds. 
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FIG. 5: Log-log plot of Binder's cumulant versus time ob- 
tained at criticality, starting from ordered initial conditions 
(mo = 1), for k=4, 5, and 6. The full lines correspond to the 
best fits obtained for t > 30 MCS. 



TABLE IL Critical exponents determined from the time de- 
pendence of the magnetization (2nd column, equation JSJ), the 
Binder cumulant (3rd column, equation(|7|l) and the logarith- 
mic derivative of magnetization (4th column, equation 0)- 
Also, * corresponds to a crude estimation performed by using 
equation Data obtained by starting the simulations from 
an ordered initial state and for the k=4, 5, and 6 segmentation 
steps of the fractal. 

i-^z z uz ^ 

4 0.0145(1) 0.556(3) 0.209(2)* 0.069(1) 

5 0.0116(1) 0.54(1) 0.167(8) 0.069(3) 

6 0.0110(1) 0.54(3) 0.162(8) 0.068(3) 



From these results the order parameter critical expo- 
nent (5 can be obtained, which appears to be less sensitive 
to the change of the segmentation step than other expo- 
nents shown in Tabled It is worth mentioning that our 
best estimation, given by /3 = 0.068(3), is significantly 
smaller than the values corresponding to the Ising model 
in d = 2 (/3 = 0.125(exact)) and in d^ ^ 1.8927 on 
SC(3,1) (/3 = 0.121(5)). 



B. Dynamic Evolution from the Disordered State 

The dynamic evolution after quenching the system to 
Tc when the simulations are started from the disordered 
state presents a weak dependence on temperature. This 



FIG. 6: Log-log plots of magnetization versus time obtained 
at critical temperature, starting from disordered initial con- 
ditions slightly modified to obtain different values of the ini- 
tial magnetization mo. Data corresponding to fc = 4(ns = 
50000), 5(71^ = 5000) and 6(n3 = 500), and different values of 
mo, which are also indicated in the figure. The inset shows the 
dependence of Q on the initial magnetization mo that allowed 
us to extrapolate the exponent S(mo 0) for each k. 



shortcoming hinders an independent estimation of Tc 
based on these measurements. However, by using the 
value of Tc obtained by means of simulations started from 
the ordered state, it is possible to obtain an independent 
evaluation of the critical exponents. 

FigureElshows the initial increase in magnetization ob- 
served for the segmentation steps 4, 5, and 6, obtained for 
different values of the small initial magnetization (mo = 
0.02, 0.04, and 0.06). The data exhibit a weak depen- 
dence on the segmentation step as can be deduced from 
the overlapping of the curves. Within the time regime 
considered, the magnetization always increases and the 
data can be fitted to a power law, as expected from equa- 
tion ||2Jl. Nevertheless, a soft curvature of the data can 
be observed for larger times due to the fact that toq is fi- 
nite and the power law is actually expected to hold in the 
mo — > limit. So, in order to determine the critical expo- 
nent we performed a fit of the data within the time inter- 
val 30-lOOMCS. As can be observed in the inset of Figure 
|S1 the exponents show a weak dependence on mo. Then, 
the exponent Q was evaluated by a linear extrapolation 
to mo — 0. So, according to our results, listed in Table 
Iml the exponent Q = 0.181(2) for the SC(4,2) fractal ap- 
pears to be the same as that for SC(3,1) {B ~ 0.1815(6)) 
and slightly smaller than the value for the Ising model in 
d = 2, given by 6* = 0.191(3) |2H. 
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FIG. 7: Log-log plots of the second moment of magnetiza- 
tion versus time, obtained at the critical temperature, start- 
ing from disordered initial conditions with mo = 0. Data 
obtained by taking fc = 4, 5, and 6. The full lines indicate the 
fit of the data according to equation ^ for t > 30. 



Figure [7| shows the time evolution of the second mo- 
ment of the magnetization and the corresponding fit ob- 
tained, within the interval 30-10'' MCS, according to 
equation |0J. The obtained values are also listed in 
Table lllll which are within the error bars independent 
of the segmentation step. The exponent 7 of the sus- 
ceptibility can be estimated assuming that hyperscaling 
law deff = ^ + ^ holds, and by combining the results 
corresponding to different initial conditions, i.e. using 
the exponents of the second moment and of the loga- 
rithmic derivative of the magnetization. The obtained 
values (see Table llllll show a convergent increase with 
k, giving 7 = 3.3(2) for fc = 6. This value is signif- 
icantly larger than those obtained for the Ising model 
in d = 2 and in SC(3,1), namely 7 — 1.75 (exact) and 
7(fc = 6) = 2.22(2), respectively. However, it is consid- 
erable smaller than the values reported from the simu- 
lations performed in equilibrium and obtained by using 
finite size scaling for the same segmentation step, which 
are 7 > 5.39[1|^. 

On the other hand, the decay of the autocorrelation 
function (see Figure 8) slightly depends on the segmen- 
tation step. The fits were carried out after disregarding 
an initial time interval oi train — 30MCS, since after that 
the power-law behavior expected from equation is ob- 
served. The obtained exponents A = def f / z — 9 are also 
reported in Table IIIII It is worth mentioning that by 
inserting the value of 6 already determined from the ini- 
tial increase in the magnetization (see Figure in the 



FIG. 8: Log-log plots of autocorrelation versus time, obtained 
at critical temperature, starting from disordered initial con- 
ditions with mo = and by taking = 4, 5 and 6. The full 
lines indicate the fit of the data according to equation ||1J for 
t > 30. 



exponent of the autocorrelation function, one can also 
calculate deff/z, as listed in the 4th column of Table 
lllll The obtained results are in full agreement with in- 
dependent determinations performed by fitting the time 
dependence of the Binder cumulant starting simulations 
from ordered states, see Table Hll 

Also, by inserting the exponent d^ff/z in the expres- 
sion of the exponent of the second moment of the mag- 
netization, given by d^ff/z — 2(3 /vz, one can obtain an 
additional estimation of (3 jvz^ as listed in the 5th column 
of Table IIIII The values are smaller than those obtained 
for the magnetization decay corresponding to the evolu- 
tion from the ordered state (see the 2nd column of Table 
nil). We attribute this difference to the small value of 



2/3 



In this 



(ijvz relative to both deff/z and {^-^ 

way, a small relative error for these exponents induces a 
large error in (3/vz^ as shown in Table ITTll 



C. Discrete Scale Invariance in Space and Time 

Discrete Scale Invariance (DSI) has extensively been 
discussed by Sornette0|. It is a weak kind of scale in- 
variance according to which an observable 0{x), which 
is a function of a control parameter x, obeys the scaling 
law 



0{x) = fi{b)0{bx), 



(10) 



under the change x bx. Here b is no longer an ar- 
bitrary real number as in the case of continuous scale 
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TABLE III: Critical exponents determined from the dynamic 
behavior of the second moment of the magnetization (2nd 
column, see equationj^), the autocorrelation (3rd column, 
see equation I^J) and the initial increase in the magnetization 
(4th column, see equation(|5J). Data obtained by starting the 
simulations from disordered initial states and for k = 4,5, 6. 
From these exponents the values of d^ff/z, and 7 are 
estimated and are listed in columns 5 and 6 respectively. 



k 



cff _ 2/3 



/3 
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4 0.532(5) 0.389(8) 0.175(2) 0.564.(8) 0.016(9) 2.55(3) 

5 0.531(4) 0.360(8) 0.182(2) 0.542(8) 0.006(9) 3.2(2) 

6 0.529(5) 0.360(9) 0.181(2) 0.541(9) 0.006(10) 3.3(2) 



invariance but it can only take specific discrete values, 
which in general have the form: 5„ — (61)", where bi is 
the fundamental scaling ratio. 

It is easy to show that if an observable 0{x) satisfies 
equation Hll)|l for an arbitrary b it necessary has to obey 
a power law of the type 0{x) = Cx"', where a is an 
exponent. But in the case of DSI the solution of equation 
(|10(l in general has the form 



0{x) = x"F 



log(fci) 



(11) 



where is a periodic function of period one. 

To illustrate that the SC(4,2) leads up to DSI, let us 
calculate the dependence of the mass for this fractal on 
the distance R to an arbitrary position Tq given by 



m(i?, ro)= / / p{r - ro)rd0dr; p{r) 5(r - n) 



N 

E 

1=1 

where i runs over all the sites in the fractal and ri is the 
position of the i-th site. Note that, in principle, the mass 
depends on the point ro we are choosing as coordinate 
origin. Figure |^ shows a log- log plot of the dependence 
of m{R, ro) on R, where ro is taken as one of the cor- 
ners of the fractal. As can be seen in the inset, the mass 
oscillating component adjusts quite well to the behavior 
expected from equation ((TT1) . The value for bi = 4.01(2) 
obtained by the fit with four terms in the Fourier ex- 
pansion coincides with the linear size of the generating 
cell, which of course is the fundamental scaling ratio. By 
choosing different positions ro, the dependence of the 
mass on R looks different than the one shown in Fig- 
ure El However, when we compute •m{R), the averaged 
mass function over all the sites ro in the lattice, we have 
find that it is still possible to factorize a periodic part 
F(log(i?)/ log(6i)) with 61 K, 4. Therefore, we conclude 
that the mass, an exclusively topological property of the 
fractal, presents a logarithmic oscillating behavior, this 
a being signature of the DSI of the fractal. 

We would like to remark that the fitted value of the 
fundamental scaling ratio bi = 4.01(2) for fc = 5 ob- 
tained from Figure © is in excellent agreement with the 
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FIG. 9: Log-log plot of the mass m{R, ro) as a function of 
R, for the segmetation step = 5 of the SC(4,2). The origin 
To has been chosen at one of the fractal corners. The inset 
shows the logarithmic oscillatory part, P(log(_R)/ log(6i)) = 
m{R, ro) / R''"), and a fit with the first four terms of the 
Fourier expansion. For the fundamental scaling ratio we have 
obtained 61 = 4.01(2). 



exact value given by bi = 4. This result suggests a fast 
convergency of the logarithmic periodic behavior towards 
the value corresponding to the mathematical fractal and 
gives us strong confidence in the values of the character- 
istic time evaluated by fitting the logarithmic periodic 
dynamics for fc = 6 (see Table (QJ). 

The observed logarithmic periodic oscillation in the de- 
cay of the magnetization (Figures |2Jl and Q) leads us 
to abandon the standard power-law decay of the form 
M(t) (X t~~ (equation ((SJ) and propose a new Ansatz 
based on equation namely the typical power law 

modulated by a logarithmic periodic function (see equa- 
tion l|SJ)). Of course, by means of this assumption we im- 
plicitly recognize that the DSI intrinsic of the susbtrate 
is somewhat capable of infiuencing the dynamic behavior 
of the overlaying ferromagnet yielding to a DSI in the 
time scale. Following this line of reasoning one has to 
admit that the fundamental scaling ratio between lengh 
scales (61) of the fractal has to translate into an intrin- 
sic ratio between time scales for the oscillations given by 
the period P. It is well known that for a critical system 
evolving towards equilibrium, characteristic spatial and 
temporal scales are linked through the dynamic exponent 
z, i.e the development of the correlations up to a length 
^ scales as ^ oc t^/^ . Therefore, going one step further 
in our speculation, we propose that the scaling ratios are 
also linked according to 



b, = pi 



(13) 
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In Appendix A we demostrate under which assump- 
tions this relation holds, but let us now show that equa- 
tion is consistent with our measurements. In fact, 
equation (|13() implies z = j^^^j^- Then, from Table we 
take log{P) — 2.16 yielding z « 3.6. Also, the logarith- 
mic derivative of the magnetization yields log{P) — 2, 
so one has z « 3.3. These rough estimations of the dy- 
namic exponent can be compared to the accepted value 
for the 2D-Ising model, namely z=2.165 plj and our pre- 
vious estimation of z = 2.55[l3 for the SC(3,1). So, all 
these results are consistent with an increasing trend of 
z when the fractal dimension of the substrate decreases, 
anticipating the occurrence of a dramatic increase in the 
slowing down, characteristic of the second-order phase 
transitions, in low dimensionality. Furthermore, inde- 
pendent measurements (see Tables ^ and IIII|I yield to 
^^^^ w 0.54, so one has de// ~ 1.9 (z « 3.6) and 
deff ~ 1.8 (z « 3.3), i.e. two figures in reasonable agree- 
ment with the fractal dimension of the SC(4,2) given by 
dn — 1.7925. On the other hand, based on both the mea- 
sured exponent ^ = 0.162(8) and the crude estimation 
of z, we could give a very rough estimation of the corre- 
lation length exponent, which should be of the order of 
J/ « 1.8, also in agreement with the trend shown for the 
2D-Ising model {ly = 1) and SC(3,1) {f « 1.39). 

IV. CONCLUSIONS 

We have studied the nonequilibrium critical dynamics 
in the short-time regime of the Ising ferromagnet em- 
bedded in a fractal substrate, namely the SC(4,2). The 
influence of the topology of the fractal on the dynamic 
evolution of some physical observables is clearly identi- 
fied through at least two effects: i) The dependence of 
the short-time regime on the segmentation step; and ii) 
the occurrence of logarithmic periodic oscillations super- 
imposed to the power-law behavior observed in both the 
decay of the magnetization and its logarithmic deriva- 
tive with respect to the reduced temperature, when the 
system is annealed from T = to Tc. We propose that 
these oscillations are a consequence of the DSI of the 
fractal substrate. 

In order to describe the decay of the magnetization 
upon annealing to we proposed an Ansatz that in- 
volves the standard power-law behavior but now modu- 
lated by a logarithmic periodic oscillatory function (see 
equation ©). This kind of function has also been used 
to describe the behavior of systems that exhibit a DSI 
in its dynamic behavior[l2[l3. By fitting the data with 
equation |(SJ we determine = 1.10(1), a figure that is 
in relatively good agreement with the values obtained by 
other authors by performing equilibrium measurements 
and by applying a finite-size scaling aT)proach|3.l20j|. 

The exponent 9 = 1.81(2) of the initial increase in 
the magnetization, determined for the segmentation step 
k = 6, appears to be the same as that for the SC(3,1) 
{6 = 0.1815(6)) but it is only slightly smaller than the ac- 



cepted value for the Ising model in d = 2 {9 = 0.191(3)), 
suggesting that this exponent is not significantly affected 
by the dimensionality of the substrate. It is worth men- 
tioning that 9 is related to xq, such that xq — z9 + — and 
that the former exponent sets the time scale for the initial 

increase in the magnetization through t„iic — TOq " ■ 

The second moment of the magnetization and the au- 
tocorrelation function obtained by starting the dynamic 
evolution from the disordered initial state follow the ex- 
pected power-law behavior. This fact allows us to per- 
form an independent determination of the exponents 
^^^j^ and which turn out to be self-consistent with 
the values obtained from simulations starting from the 
ordered state. 

We would like to remark that the critical exponents 
corresponding to all observables measured show a con- 
vergent trend when the segmentation step is increased, 
indicating that the largest segmentation used in this work 
(fc = 6) would be considered a good approximation of the 
mathematical fractal. 

Our estimation (3 = 0.068(3) indicates that the order 
parameter critical exponent for the SC(4,2) is smaller 
than that corresponding to the Ising model in d = 2 (/? = 
1/8), which is almost the same as that reported for the 
case as the SC(3,1) {(3 = 0.121(9)). On the other hand, 
our estimation for the exponent of the susceptibility (7 = 
3.3(2)) is significantly larger than those corresponding to 
the 2D-Ising magnet (7 — 1.75)) and the SC(3,1) for the 
same k (7 = 2.22(2)). Nevertheless, it is notably smaller 
than the values reported from e quili brium simulations 
using finite-size scaling (> 5.39) [lll^. 

We hope that the present work represents an extensive 
attempt to not only numerically characterize the rele- 
vant critical properties of the Ising model on a fractal 
substrate, but also to give clear evidence that the DSI of 
the underlying fractal structure infiuences the dynamic 
evolution of an Ising magnet. Of course, we have dis- 
cussed the difficulties one encounters when dealing with 
these systems, but even more importantly our study al- 
lows us to clearly identify and characterize the influence 
of the substrate topology on the power-law behavior ex- 
pected for the nonequilibrium critical dynamics in the 
short-time regime. The DSI of the fractal shows up dra- 
matically in the dynamic behavior of the magnetization 
decay by causing the occurrence of logarithmic periodic 
oscillations. 



V. APPENDIX A 

In this appendix we show that equation H13|) can be 
straightly proven for a certain observable M, if one as- 
sumes: 

(i) M{t) obeys a time DSI of the form: 

M{t) ^ tx{P„)M{Pr,t), Pn^P" (14) 

where P is the period and n is an integer. 



(ii) The correlation lenght scales as 

eoc<l/^ (15) 

(iii) The same observable M, as a function of the cor- 
relation lenght, obeys a spatial DSI of the form: 

M{0 = Kbn)M{bnO, K^h^l (16) 

where bi is the fundamental scaling ratio between 
lenght scales and M is using to denote: M{S) — 

Mm]. 

Assumption [ii) implies t — a^^ for some constant a. 
Then by replacing on the right-hand side of equation 114(1 
one obtains 

M{t) = tJ,{P'')M{P''a^'^) = /x(P")M[a(P^O^], (17) 
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and writing M{t) as M{^) this equality becomes 

Alio = M(P")M(P^e), (18) 

The comparison with equation (|16|l leads us to the the 
following equalities 



and therefore, = hi. 
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